Generalization of the Beck-Cohen superstatistics 
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Generalized superstatistics, i.e., a "statistics of superstatistics," is proposed. A generalized su- 
perstatistical system comprises a set of superstatistical subsystems and represents a generalized 
hyperensemble. There exists a random control parameter that determines both the density of en- 
ergy states and the distribution of the intensive parameter for each superstatistical subsystem, 
thereby forming the third, upper level of dynamics. Generalized superstatistics can be used for non- 
stationary nonequilibrium systems. The system in which a supercritical multitype age-dependent 
branching process takes place is an example of a nonstationary generalized superstatistical system. 
The theory is applied to pair production in a neutron star magnetosphere. 
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I. INTRODUCTION 

A large variety of complex nonequilibrium systems ex- 
hibit spatiotemporally inhomogeneous dynamics. Such 
systems are often characterized by hierarchical structures 
of dynamics. The hierarchy can be formed by the decom- 
position of the system dynamics into different dynam- 
ics on different spatiotemporal scales, which are largely 
separated from each other. In this case, the statistical 
properties of the system can be effectively described by a 
superposition of several statistics, or a "superstatistics." 

Superstatistics has been formulated in Ref. [J to con- 
sider nonequilibrium systems with a stationary state and 
intensive parameter fluctuations. Some preliminary con- 
cepts have been anticipated earlier in Refs. ^2i-t^]. Super- 
statistical systems are characterized by the existence of 
an intensive parameter j3 that fluctuates on a much larger 
time scale than the typical relaxation time of the local 
dynamics. If a given system can be thought of as a col- 
lection of many small spatial cells, then the inverse tem- 
perature in a cell is typically taken as such a parameter. 
However, more general interpretations of the intensive 
parameter are possible. Sufficient time scale separation 
between two relevant dynamics within the complex sys- 
tem allows one to qualify superstatistics as a form of slow 
modulation |5j]. 

Superstatistics is applicable to various complex sys- 
tems. Its applications include, among others, cosmic- 
ray energy spectra and electron-positron pair annihila- 
tion 01, the world line representations of J^eynman 
propagators for spin-0 and spin- 1/2 particles |8|, an ex- 
tension of the random matrix theory covering systems 
with mixed regular-chaotic dynamics [9l-[lll|. nonstation- 
ary dynamical processes with time- varying multiplicative 
noise exponents [l3| , Markovian systems without detailed 
balance [l^, a mesoscopic approach to the problem of 
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Brownian motion |14|, m odels of the metastatic cascade 
in cancerous systems [l^ , complex network s IIGII , ecosys- 
tems driven by hydroclimatic fluctuations [17|, pattern- 
forming systems [l8| . solar flares [l^, share price fluc- 
tuations |20l - [23j , the statistics of train departure delays 
[2^ , wind velocity fluctuations [1^ , and many interesting 
applications in hydrodynamic turbulence [23, l26l - [30j . 

Let us suppose for a while that we have a complex 
nonequilibrium system described by superstatistics. One 
of the main problems is the determination of the dis- 
tribution /(/3) of the intensive parameter /3. It can be 
found using the maximum entropy principle [3l| . More- 
over, /(/3) can be considered as a function of a set of 
additional control parameters {A^}, / = f{j3,{\i\) [3^. 
These parameters emerge as Lagrange multipliers deter- 
mined by maximizing the Boltzmann-Gibbs-Shannon en- 
tropy of /(/3, {Xi}) under certain constraints. Let us re- 
vert to the parameter /?. In each cell of the supersta- 
tistical system, we have the Gibbs canonical distribution 
Pg[E\I3) = e-^^/Z(/3), where Z(/3) is the partition func- 
tion. This distribution can be obtained by maximizing 
the Boltzmann-Gibbs-Shannon entropy associated with a 
normalized distribution p(i?|/3), where /3 is the Lagrange 
multiplier corresponding to the mean energy constraint. 
We see that, in general, both (3 and {A^} have a similar 
nature. However, in superstatistics, /3 fluctuates, but the 
control parameters are constant. 

We can pose the following question: What statistics 
appears in the case of fluctuating control parameters? 
The aim of this paper is to develop the generalization of 
the Beck-Cohen superstatistics that allows one to prop- 
erly consider such fluctuations. These form a separate, 
upper level of dynamics, thus leading to multiscale super- 
statistics. The possible existence of such superstatistics 
has been presumed in Ref. (3T| . 

This paper is organized as follows. In Sec. II general- 
ized superstatistics is developed. In Sec. Ill supercritical 
multitype age-dependent branching processes are consid- 
ered as an example of generalized superstatistical sys- 
tems. In Sec. IV generalized superstatistics is applied to 
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electron-positron pair production in a neutron star mag- 
netosphere and the particle energy distribution is calcu- 
lated. In Sec. V the main results are summarized. 



II. GENERALIZED SUPERSTATISTICS 

Let us consider a nonequilibrium system that consists 
of a set of nonequilibrium subsystems described by their 
superstatistics. This system can be referred to as a "gen- 
eralized superstatistical system." Let us define a ran- 
dom variable ^ that determines the properties of each 
nonequilibrium superstatistical subsystem. The role of ^ 
in determining the properties of the subsystems is simi- 
lar to that of /3 in determining the properties of the cells 
that constitute a superstatistical subsystem. In other 
words, ^ is a control parameter, which determines the 
form of the superstatistical distribution for each subsys- 
tem. Here, we do not restrict ourselves to considering ^ 
as a scalar random variable. In principle, ^ is allowed to 
be a multidimensional random vector. The distribution 
of ^ is characterized by a probability density c(^), which 
is normalized, / c{£_)d^ — 1. 

We can look at the generalized superstatistical system 
using the hyperensemble approach [33]. Developed in 
the context of nonequilibrium systems, this approach is 
consistent with the standard theory of equilibrium sta- 
tistical physics [si^. Superstatistics can be considered as 
the theory of hyperensembles [35|. Each superstatisti- 
cal subsystem of the system represents a hyperensemble, 
i.e., a mixture of canonical ensembles. In turn, the sys- 
tem as a whole can be described by a mixture of the hy- 
perensembles corresponding to the subsystems. In other 
words, the generalized superstatistical system represents 
a "generalized hyperensemble," i.e., an ensemble of hy- 
perensembles. 

There are three levels of dynamics in the system under 
consideration: The first level is the level of fast dynamics 
in a cell. The second level is the level of superstatistical 
dynamics in a subsystem. The existence of an additional 
dynamics described by ^ forms the third level in the hier- 
archy of dynamics and allows us to refer to the statistics 
of the generalized superstatistical system as a "general- 
ized superstatistics," i.e., a statistics of superstatistics. 

Consider a superstatistical subsystem of the general- 
ized superstatistical system. Let r(£'|^) be a nondecreas- 
ing function representing the number of states with en- 
ergy less than E. The random variable f determines the 
density of energy states for the subsystem. 



for each cell of the subsystem is 



dT{E\0 
dE ' 



The inclusion of the density of states in the superstatisti- 
cal description is significant [s^. When writing dT{E\£^) 
in integrals, we will imply integration over E so that 
fir(£;|^) = g{E\^)dE. The Gibbs canonical distribution 
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where 



(1) 



is the partition function. By analogy with the approach 
used in Ref. [3^ , we do not include the density of states 
in the definition of the energy distribution ([Ij ; therefore, 
the normalization condition is / pG{E\f3,£,)dr{E\£^) = 1. 
We also assume that ^ determines the distribution /(/3|^) 
of the intensive parameter /3. This distribution is normal- 
ized, //(/3|^)d/3 = 1. Note that g{E\^) and are 
not necessarily statistically dependent variables. These 
will change independently if we choose ^ = (^1,^2), with 
^1 and ^2 being independent random variables, and set 

g{E\0 ^ 9{E\^i) and /(/3|0 = /(/3|6)- 

The superstatistical distribution for the subsystem is 
given by 



(2) 



with the normalization condition being 
/ p{E\£^)dT{E\£^) = 1. If ^ is a nonrandom vector, 
then this distribution reduces to the ordinary supersta- 
tistical distribution f]\. By averaging Eq. © over the 
fluctuating ^, we immediately obtain the generalized 
superstatistical distribution 

aiE) = I piE\OgiE\Ocim, (3) 

which is normalized, J a{E)dE ~ 1. 

III. BRANCHING PROCESSES 

Now the following question arises: Do there exist any 
nonequilibrium systems, whether stationary or not, that 
can be described by generalized superstatistics? 

Let us consider a many-particle system composed of 
particles of n types. Each type-? particle, which will also 
be denoted by T^, 1 ^ i ^ n, has a random lifetime 
with a probability distribution function Gi(r). In other 
words, Gilr) is the probability that the lifetime of a given 
type-i particle does not exceed r. At the end of its life 
the particle decays into a random number of particles of 
several types. Specifically, at the moment of its decay 
the particle produces LUj ^ type-j particles of age zero, 
1 f^j f^n: 



(4) 



Thus we have a multitype age-dependent branching pro- 
cess, the so-called multitype Sevast'yanov process [371] . 
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The transformation ^ is described by the generating 
function 



/i,(t,x)=^P-(tK^...<" 



where p"(t) is the conditional probability of the trans- 
formation (|4]) given that the particle decays at age r, 
Lo = (wi, . . . jWn) is an n-dimensional vector with non- 
negative integer components Wj, and x = (xi, . . . ,Xn) 
is an n-dimensional vector with complex components Xj 



such that 



^ 1. If Xi 



1, we will write 



X = 1. The normalization condition reads /ij(r, 1)=1. 

The mean number of type-j particles that appear upon 
the decay of a type-« particle is given by 



Ai-j — 



where 



aij{T)dGi{T) 



dhi{T,x) 



(5) 



dxi 



is the same mean given that the particle decays at age t. 
We will assume that the n x n matrix A ~ \\Aij\\ with 
components ^ Aij < oo is irreducible, or indecompos- 
able, i.e., the index set {1, . . . , n} cannot be divided into 
two disjoint nonempty sets 5*1 and 5*2 such that Aij — 
for all i £ Si and all j € S2 (see, e.g., Ref. [SS] for de- 
tails). Moreover, we will also assume that the Perron 
root of A, i.e., the maximum positive real eigenvalue of 
A, is greater than one. Thus we deal with the indecom- 
posable supercritical branching process [s^. Physically, 
this means that, first, a particle of a given type poten- 
tially has descendants, either direct or distant, of any 
type and, second, the number of particles in the system, 
on average, progressively increases. 

We see that the whole system is composed of n sub- 
systems, the ith subsystem comprising type-i particles. 
Obviously, the subsystems interact with each other in the 
sense that the decay of a particle in one subsystem leads 
to the creation of particles in other subsystems. The 
number of particles both in the whole system and in the 
subsystems is not constant and changes with time. Thus 
we have a nonstationary nonequilibrium situation. 

Let fi{t) = [fii{t), . . . , fJ-n{t)] be an n-dimensional vec- 
tor with non-negative integer components. Each fj.j{t) 
is a random variable that yields the number of parti- 
cles in the jth subsystem at time t. The age-dependent 
branching process is characterized by a set of generating 
functions 



(6) 



where Pi{t) is the conditional probability that — uj 
given that there is one particle in the ith subsystem and 
none in the other subsystems at time zero. The generat- 
ing functions Fi{t,x) 1 ^ i ^ satisfy a system of 



nonlinear integral equations ■ Essential for us here is 
the probabilistic nature of the number of particles both 
in the whole system and in each subsystem. In general, 
we cannot predict a priori what the number of particles 
is in each subsystem at a given time, but instead we can 
focus on the average behavior of the generalized hyper- 
ensemble. 

The mean number of particles in the jth subsystem at 
time t given that the branching process has started with 
one type-? particle at time zero is 



A,,{t)^ 



dF,{t,x) 



dxi 



To analyze the long-run behavior of these means, let us 
define the Laplace-Stieltjes transforms 



where 



Lij (a) 



e-"^dG,,(r), 



aij{T)dGi (t). 



(7) 



Let A (a) be the Perron root of the n x n matrix L{a) — 
\\Lij(a)\\. Choose an a such that X{a) ~ 1 and define an 
n-dimensional left eigenvector v — (ui,...,u„) of L{a) 
with positive real components Vj that satisfy 



E 



ViLij{a). 



Then the asymptotic behavior of the mean numbers of 
particles is expressed as 



A,,{t) 



t 



00, 



where 



e-°'^[l-Gj{T)]dT 



(8) 



(9) 



and Ci is a positive constant [40|. Interestingly, Aij{t) 
behaves in a similar way for each subsystem. Since we 
consider the supercritical case, a always exists and is pos- 
itive. Therefore, the mean number of particles in each 
subsystem increases exponentially. Note that the exis- 
tence of a positive a is equivalent to the condition that 
the Perron root of A is greater than one [s^, . 

In the long-term run, the limiting probability that a 
given particle belongs to the ith subsystem can be written 
as 



(10) 



En 



Note that tt^ is independent of the type of the primary 
particle with which the branching process has started. 
Moreover, nonstationary though the situation is, the lim- 
iting probability is stationary. 
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The hmiting age distribution for the ith subsystem is 



L^ir) 



/;e-°"[l-G,(u)]dM 
J-e-""[l-G,(u)]du' 



(11) 



It is found by dividing the ith subsystem into two sub- 
systems, the first comprising the type-i particles that de- 
cay at age u ^ t and the second comprising the rest, 
and applying the technique described above to the new 
system. Equation (|lip yields the probability that a ran- 
domly chosen type-i particle decays at age u ^ t. In the 
case of a single type of particles, Li{T) reduces to the 
limiting age distribution for the one-dimensional super- 
critical Sevast'yanov process [4l|, which coincides with 
the classical age distribution for the one-dimensional su- 
percritical Bellman-Harris process when t he p robability 
pfir) is independent of r, p^(t) = pf [H, ig. 

Now we can calculate the particle energy distribution. 
Consider a type-i particle of age r. In general, its en- 
ergy can be considered as a random variable character- 
ized by a conditional probability density Wi{E\T). More 
specifically, the probability that the energy of a type-i 
particle of age r lies in a small interval dE around E 
is 'Wi{E\T)dTi{E), where ^i{E) is the number of energy 
states with energy less than E. The normalization con- 
dition is / Wi{E\T)dTi{E) = 1. The energy probability 
density for the ith subsystem becomes 



W^{E\T)dL,{T), 



(12) 



with the normalization condition being J pi{E)dT i{E) — 
1. 

The described system can be considered as a general- 
ized superstatistical system. The control parameter ^ is a 
discrete random variable that takes on values {1, . . . , n} 
and yields the number of the subsystem to which a ran- 
domly chosen particle belongs. In other words, ^ corre- 
sponds to the particle type. It has the discrete prob- 
ability distribution {tti, . . . , 7r„}, where tt^ is given by 
Eq. ([TU]). It remains to find the distribution /i(/3) of 
the fluctuating parameter j3 for each subsystem. Let 

/•oo 

g{s)=Z[f{x)]{s)^ / e-'-f{x)dx 
Jo 

be the Laplace transform of a function f{x), with f{x) — 
£,~^[g{s)]{x) being the corresponding inverse Laplace 
transform. The energy distribution for the ith subsys- 
tem can be expressed as 



P^{E) 



MP) 



where Zi(/3) is the partition function. Then 
MP) = Z,{P)£-'[p.,iE)m. 



IV. AN EXAMPLE: PAIR PRODUCTION IN A 
NEUTRON STAR MAGNETOSPHERE 

New nonstationary cosmic radio sources associated 
with neutron stars, viz., intermittent pulsars [4^ and 
rotating radio transients (RRATs) [i^, have been dis- 
covered recently. The essential feature of these sources is 
their long "silence," when we do not observe any radio 
emission from them. Since an electron-positron plasma 
outflowing from the magnetosphere of a neutron star is 
responsible for the observable radio emission, the plasma 
generation can be switched off for some time. In this 
case, the absorption of a high-energy photon in the in- 
ner neutron star magnetosphere triggers nonstationary 
cascade pair production [46|, which, in turn, results in 
the formation of a "lightning" [13]. Such lightnings can 
manifest themselves as radio bursts from RRATs [48| . 
In a lightning, we deal with an ultrarelativistic electron- 
positron plasma. The properties of the emission from 
electrons and positrons are determined by their energies. 
Therefore, it is important to find the energy distribution 
of particles. This can be done using generalized super- 
statistics. 

Let us characterize the energy of a charged particle, 
either an electron or a positron, by its Lorentz factor 
7(t). The particle is efficiently accelerated by a longi- 
tudinal electric field so that 7(r) eventually reaches 
a stationary value 70, which is ^10^ in a vacuum neu- 
tron star magnetosphere 49] . The electron and positron 
of each produced pair, though ultrarelativistic, initially 
have Lorentz factors much less than 70. At the initial 
stage of acceleration j(t) increases linearly with time, 
7(t) « E^^T. Here, we use a dimensionless system of units 
(see, e.g., Ref. [3). By contrast, when t approaches 
■'0 = 7o/£'||j the radiation forces come to the fore, and a 
need arises to use the Dirac-Lorent z eq uation to consider 
the particle dynamics properly p9l - [5l| . 

Thus, we can define two types of particles: A type-1 
particle can be efficiently accelerated by the electric field 
since the radiation friction is negligible. A type-2 par- 
ticle, in contrast to a type-1 particle, is not accelerated 
by the electric field because of the electrodynamic self- 
action effects and has the constant Lorentz factor 70. In 
the former case the particle does not efficiently produce 
secondary pairs, but in the latter case it does at a rate 
Q [4^. Note that the plasma generation, along with the 
accompanying radio emission, is not suppressed even in 
ultrahigh magnetar magnetic fields ^3] . The Lorentz 
factors of type-1 and type-2 particles as functions of their 
ages become 



72 (t) 



--E\\r, 
- 70, 



s=: r < To, 
r < 00. 



(13a) 
(13b) 



The particles of each produced pair, though moving 
independently of each other, can conveniently be consid- 
ered as a whole, and type-1 and type-2 pairs, which will 
be denoted by Ti and T2, respectively, are defined by 
analogy with individual particles. 
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Now we can write the following transformations: 

Ti ^ T2, 

Ta ^Ti +T2. 

The generating functions are 

hi{T,x) = X2, 
/12(t, x) = X1X2, 



where x = [xi, X2)- The 2x2 matrix a 
( 5 J ) . The lifetime distribution functions are 



aij 1 1 becomes 



Gi(r) 
G2{r) 



To), 

-2Qt 



(14a) 
(14b) 



where 9{x) is the Heaviside function. The matrix A de- 
fined by Eq. ([5|) coincides with a and is indecomposable 
since its off-diagonal elements are positive. The Perron 
root of A is equal to (1 + \/h)/2 and greater than 1. 
Therefore, we deal with the indecomposable supercriti- 
cal two-type branching process. 

Now we should find a proper value of a as described 
after Eq. ([7]). Interestingly, we can do this without di- 
rectly finding the Perron root of L{a). We have recently 
shown [46| that pair production in the system under con- 
sideration is asymptotically described by the equation 



dN{t) 
dt 



2Q''*'N{t), 



(15) 



where N{t) is the number of electron-positron pairs at 
time t, Q"^^ = /2tq is the effective pair production 
rate, and N!^l^ satisfies 

iV^f =lnAr.„-lniV;ff, 

where Nrg — 2Qto is the number of particles created by 
a fully accelerated particle in time tq. Note that N^^ can 
be equivalently expressed through the Lambert function 
W{x) (see, e.g., Ref. [H^l for its definition) as 

and hence is positive since W{x) is positive for any posi- 
tive x. From Eq. (fT5|) it follows that the total number of 
electron-positron pairs in the system increases exponen- 
tially. Comparing Eq. to Eq. (|H]) yields 



a = 2Q 



eff 



(16) 



which is positive, as it must be in the supercritical case. 

The energy of any type-1 particle is less than 70, while 
that of any type-2 particle is 70. Therefore, it is natural 
to choose the following density of states for the subsys- 
tems: 



51(7) = 1 - ^(7-70) 
52(7) = (5(7- 70), 



(17a) 
(17b) 



where S{x) is the delta function. We see from Eq. 
that the energies of particles are nonrandom functions 



of their ages, hence the conditional energy distributions 
have the form 



wi{j\t) =(5[7-7i(t)], 
■^2(7^) = 1- 



(18a) 
(18b) 



Using Eqs. (HH, dH]), dH]), and JTH]), we obtain the en- 
ergy distributions for the subsystems: 



Pi (7) 



E\i 1 - e-"-^o ' 



(19a) 
(19b) 



^2(7) = 1- 

The corresponding intensive parameter distributions are 
/2(/3) = SiP). 

Let us find the probability distribution of the random 
control parameter ^ that corresponds to the type of a 
randomly chosen particle. From Eqs. ^ and (|14p we 
have 



Wi 
W2 



1 - e^^^o 



1 



(20a) 
(20b) 



a + 2Q 

The components of the left eigenvector v — {vi, V2) of 
/ e-"^° \ 



can be chosen as follows: 

Vl = 1, 



(21a) 
(21b) 



Using Eqs. ([TH), ([20]), and we obtain the prob- 
ability TTj that a randomly chosen particle is of type ^, 
e = l,2: 



771 ~ 1 — 



2g' 



71'2 



a 

2Q' 



(22a) 
(22b) 



Note that tt2 may be interpreted as the probability that 
the particle significantly contributes to pair production. 

Finally, Eqs. ©, dH]), (HZ]), (Ull), and (gl]) allow us to 
obtain the generalized superstatistical distribution 



^(7) = 7^ ^(7 - 70) + [1 - 0{l - 7o)]^e-"^/^ii 



2Q 



(23) 



which represents the energy distribution of ultrarelativis- 
tic electrons and positrons. 

The particle energy distribution in the neutron star 
magnetosphere cannot be observed directly. Moreover, 
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it is unclear how to infer the particle energy distribu- 
tion for RRATs from currently available observational 
data. On the one hand, it is difficult to detect radio emis- 
sion from RRATs because of its sporadic nature, and one 
has to carry out single-pulse searches, which require sig- 
nificant radio telescope resources with long observation 
times [5^ [5^. On the other hand, it is unclear how the 
form of the distribution correlates with the properties of 
radio emission. Although the parameters of an electron- 
positron plasma generated in the RRAT magnetosphere 
imply the possibility of the efficient generation of radio 
emission j47l . l48j , the detailed mechanism is unknown at 
the moment. Therefore, the distribution ((23)) is a theo- 
retical prediction. 



V. CONCLUSION 

I have developed generalized superstatistics, i.e., a 
statistics of superstatistics. Based on the concept of fluc- 
tuating control parameters, it is defined for a general- 
ized superstatistical system, which consists of a set of 
nonequilibrium superstatistical subsystems. Such a sys- 
tem can be considered as a generalized hyperensemble, 
i.e., an ensemble of hyperensembles. A fluctuating con- 
trol parameter, which may be a multidimensional random 
vector, determines the form of the superstatistical distri- 
bution for each subsystem. In general, it determines not 



only the distribution of the intensive parameter but also 
the density of energy states. In a generalized supersta- 
tistical system, there appears the third, upper level in 
the hierarchy of dynamics besides two levels that exist 
in each superstatistical subsystem. Interestingly, gen- 
eralized superstatistics can be applied to nonstationary 
nonequilibrium systems. As an example of a nonstation- 
ary generalized superstatistical system, I have studied the 
system in which the supercritical multitype Sevast'yanov 
process takes place. In this system, the transformation of 
particles of several types occurs, and the number of par- 
ticles both in the whole system and in each subsystem 
increases exponentially. In addition, I have considered 
an astrophysical application of generalized superstatis- 
tics and obtained the energy distribution of ultrarela- 
tivistic electrons and positrons produced in a neutron 
star magnetosphere. I have found the probability that a 
randomly chosen particle significantly contributes to the 
production of secondary electron-positron pairs. 

It seems interesting to find general principles that 
might allow us to obtain the distribution of both the in- 
tensive parameter and the control parameter. Note that 
in the case of ordinary superstatistics, the main approach 
to finding the intensive parameter distribution is maxi- 
mizing either a generalized entropy or the Boltzmann- 
Gibbs-Shannon entropy under certain constraints [sol - 
[ssl . Is?! . I believe that the same approach might be used 
in the case of generalized superstatistics. However, this 
problem requires a separate study. 



C. Beck and E. G. D. Cohen, Physica A 322, 267 (2003). 
G. Wilk and Z. Wlodarczyk, Phys. Rev. Lett. 84, 2770 
(2000). 

C. Beck, Phys. Rev. Lett. 87, 180601 (2001). 

A. G. Bashkirov and A. D. Sukhanov, J. Exp. Theor. 

Phys. 95, 440 (2002). 

P. AUegrini, F. Barbi, P. Grigolini, and P. Paradisi, Phys. 

Rev. E 73, 046136 (2006). 

C. Beck, Physica A 331, 173 (2004). 

C. Beck, Eur. Phys. J. A 40, 267 (2009). 

P. Jizba and H. Kleinert, Phys. Rev. D 82, 085016 (2010). 

A. Y. Abul-Magd, Physica A 361, 41 (2006). 

A. Y. Abul-Magd, Phys. Rev. E 73, 056119 (2006). 

A. Y. Abul-Magd, B. Dietz, T. Friedrich, and A. Richter, 

Phys. Rev. E 77, 046202 (2008). 

S. M. Duarte Queiros, Braz. J. Phys. 38, 203 (2008). 
I. Lubashevsky, R. Friedrich, A. Heuer, and A. Ushakov, 
Physica A 388, 4535 (2009). 

R. F. Rodriguez and I. Santamaria-Holek, Physica A 385, 
456 (2007). 

L. Leon Chen and C. Beck, Physica A 387, 3162 (2008). 
S. Abe and S. Thurner, Phys. Rev. E 72, 036102 (2005). 
A. Porporato, G. Vico, and P. A. Fay, Geophys. Res. 
Lett. 33, L15402 (2006). 

K. E. Daniels, C. Beck, and E. Bodenschatz, Physica D 
193, 208 (2004). 

M. Baiesi, M. Paczuski, and A. L. Stella, Phys. Rev. 
Lett. 96, 051103 (2006). 



[20] M. Ausloos and K. Ivanova, Phys. Rev. E 68, 046122 

(2003) . 

[21] P. Jizba and H. Kleinert, Phys. Rev. E 78, 031122 (2008). 
[22] C. Anteneodo and S. M. Duarte Queiros, J. Stat. Mech.: 

Theory Exp. 2009, P10023. 
[23] E. Van der Straeten and C. Beck, Phys. Rev. E 80, 

036108 (2009). 
[24] K. Briggs and C. Beck, Physica A 378, 498 (2007). 
[25] S. Rizzo and A. Rapisarda, AIP Conf. Proc. 742, 176 

(2004) . 

[26] A. M. Reynolds, Phys. Rev. Lett. 91, 084503 (2003). 
[27] S. Jung and H. L. Swinney, Phys. Rev. E 72, 026304 

(2005) . 

[28] C. Beck, E. G. D. Cohen, and H. L. Swinney, Phys. Rev. 

E 72, 056133 (2005). 
[29] C. Beck, Phys. Rev. Lett. 98, 064502 (2007). 
[30] S. Abe, Phys. Rev. E 82, 011131 (2010). 
[31] S. Abe, C. Beck, and E. G. D. Cohen, Phys. Rev. E 76, 

031102 (2007). 

[32] E. Van der Straeten and C. Beck, Phys. Rev. E 78, 

051101 (2008). 
[33] G. E. Crooks, Phys. Rev. E 75, 041119 (2007). 
[34] J. Naudts, AIP Conf. Proc. 965, 84 (2007). 
[35] S. Abe, Cent. Eur. J. Phys. 7, 401 (2009). 
[36] F. Sattin, Eur. Phys. J. B 49, 219 (2006). 
[37] B. A. Sevast'yanov, Theory Probab. Appl. 9, 521 (1964). 
[38] F. R. Gantmacher, Applications of the Theory of Matrices 

(Interscience, New York, 1959). 



7 



[39] V. A. Vatutin and A. M. Zubkov, J. Sov. Math. 39, 2431 

(1987). 

[40] B. A. Scvast'yanov and V. P. Chistyakov, Theory Probab. 

Appl. 16, 199 (1971). 
[41] A. Yakovlev and N. Yanev, Stat. Probab. Lett. 77, 503 

(2007). 

[42] K. B. Athreya and N. Kaplan, Ann. Probab. 4, 38 (1976). 

[43] T. Kuczek, Ann. Probab. 10, 252 (1982). 

[44] M. Kramer, A. G. Lyne, J. T. O'Brien, C. A. Jordan, 
and D. R. Lorimcr, Science 312, 549 (2006). 

[45] M. A. McLaughlin, A. G. Lyne, D. R. Lorimer, 
M. Kramer, A. J. Faulkner, R. N. Manchester, J. M. 
Cordes, F. Camilo, A. Possenti, I. H. Stairs, G. Hobbs, 
N. D'Amico, M. Burgay, and J. T. O'Brien, Nature (Lon- 
don) 439, 817 (2006). 

[46] Ya. N. Istomin and D. N. Sob'yanin, J. Exp. Theor. Phys. 
113, 592 (2011). 

[47] Ya. N. Istomin and D. N. Sob'yanin, J. Exp. Theor. Phys. 
113, 605 (2011). 

[48] Ya. N. Istomin and D. N. Sob'yanin, Astron. Lett. 37, 
468 (2011). 



[49] Ya. N. Istomin and D. N. Sob'yanin, J. Exp. Theor. Phys. 

109, 393 (2009). 
[50] Ya. N. Istomin and D. N. Sob'yanin, Astron. Rep. 54, 

338 (2010). 

[51] Ya. N. Istomin and D. N. Sob'yanin, Astron. Rep. 54, 
355 (2010). 

[52] Ya. N. Istomin and D. N. Sob'yanin, Astron. Lett. 33, 
660 (2007). 

[53] Ya. N. Istomin and D. N. Sob'yanin, AIP Conf. Proc. 

983, 298 (2008). 
[54] S. Asgarani and B. Mirza, Physica A 387, 6277 (2008). 
[55] M. A. McLaughlin, A. G. Lyne, E. F. Keane, M. Kramer, 

J. J. Miller, D. R. Lorimer, R. N. Manchester, F. Camilo, 

and I. H. Stairs, Mon. Not. R. Astron. Soc. 400, 1431 

(2009). 

[56] E. F. Keane, M. Kramer, A. G. Lyne, B. W. Stappers, 
and M. A. McLaughlin, Mon. Not. R. Astron. Soc. 415, 
3065 (2011). 

[57] C. TsalUs and A. M. C. Souza, Phys. Rev. E 67, 026106 
(2003). 



